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Abstract. It has been recently reported that a granular mixture in which grains 

■ differ in their restitution coefficients presents segregation: the more inelastic par- 
ticles sink to the bottom. When other segregation mechanisms as buoyancy and 

, , the Brazil nut effect are present, the inelasticity induced segregation can com- 

^"b^ ■ pete with them. First, a detailed analysis, based on numerical simulations of two 

^ ' dimensional systems, of the competition between buoyancy and the inelasticity 

induced segregation is presented, finding that there is a transition line in the 
. parameter space that determines which mechanism is dominant. In the case of 

neutrally buoyant particles having different sizes the inelasticity induced segre- 
gation can compete with the Brazil nut effect (BNE). Reverse Brazil nut effect 

■ (RBNE) could be obtained at large inelasticities of the intruder. At intermediate 
^ ' values, BNE and RBNE coexist and large inelastic particles are found both near 
O the bottom and at the top of the system. 

1 Introduction 

It is known from decades that when a mixture of two types of granular particles are externally 
, excited by means, e.g. a vibration, grains of different size, shape, mass or mechanical properties 

OO ' can mix or segregate. For instance in a mixture of small grains and one large intruder, when 

I vertically vibrated, the intruder can go up [T] or down |51 -the so called direct and reverse Brazil 

. nut effects, respectively. This effect hat has been studied in many papers (see, e.g. Ref. [3] and 

references therein) , where different mechanisms have been proposed to explain the segregation 
phenomenon. For example, percolation, arching, void filling, convection, interstitial fluid effects, 
condensation, and some more have been invoked to be the cause of segregation. All these effects 
are explained in Ref. [51, although the problem is not completely understood yet. When both 
• rH , species have similar sizes (but possibly different) we can select few cases where a variety of seg- 

/\ ' regation mechanisms and scenario appear [HOIB]. For instance, in Ref. [7] particles of different 

^ , masses, radii and restitution coefficients are placed in a dish which is horizontally vibrated, 

finding complete segregation. Segregation is also found in the same geometry when the grains 
have different friction coefficient with the base [5]. Under horizontal swirling, radial segrega- 
tion of particles of different sizes has been observed 0. In avalanches, grains of different shape 
segregate in stripes [10|; in partially filled rotating drums, axial size segregation develops [Tl] . 
In two dimensional systems under gravity, sinusoidally vibrated, clustering has been observed 
|12) . This segregation effect can be modulated by using non-sinusoidal vibration . 

In some of the cases mentioned above the grain species differ on the friction or restitution 
coefficient. However few papers have studied segregation when this is the only difference between 
grains. One of these cases is Ref. [13], where a mixture of spheres that only differ in friction 
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coefficients (static, dynamic and rolling) is horizontally vibrated. They find complete mixing 
-that is, no segregation- for a fiat plate while segregation is only observed when the plate was 
slightly inclined. Therefore, these results contradict the previously mentioned ones. 

In a theoretical approach Ref. [15] constructs the hydrodynamic equations from the Boltz- 
mann equation, finding inelasticity induced segregation. The authors explain the phenomenon 
as a consequence of the temperature gradient in the system induced by inelastic collisions, and 
relate the concentration gradient with the temperature gradient. In the same spirit, Ref. |16j 
studies the low density hydrodynamics of a mixture in the so called tracer limit, i.e. where the 
concentration of one of the components tends to zero. Among other results, they find that the 
temperature ratio of both species must be a constant. This constant value had been already 
measured by two experimental groups, both in 3D |17) and 2D |18] and by means of computer 
simulations in 3D [T^ and 2D [5D]. They found that the temperature ratio reaches a constant 
value in the regions of the system where the density is low. No tracer limit is needed. Gener- 
alization to moderate density has been done by Garzo [H] based on a kinetic approach using 
the Enskog equation for dissipative hard spheres. 

Recently, it has been shown numerically that in moderately dense vibrofiuidized granular 
matter, inelasticity induced segregation takes place [32] . When a mixture of grains of equal 
size and mass, but differing in their restitution coefficient in put in a vibrating box, the more 
inelastic grains sink, segregating partially with the other species. In this paper we study if 
this inelasticity induced segregation can compete with two other known mechanisms of grain 
segregation: buoyant forces when grain densities differ and the Brazil nut effect when their sizes 
differ. 

The structure of this paper is as follows. In Section [5] we describe the system under con- 
sideration. Section |3] reproduces the results presented in Ref. [22 regarding the macroscopic 
segregation of two species differing only in their restitution coefficient. Section |4] analyzes the 
competition with buoyant forces, when considering grains of equal size, but different masses 
and restitution coefficients. Section [5] studies the competition with the Brazil nut effect when 
grains of the two species have the same density, but differ in size and restitution coefficients. 
We conclude with Section |5] summarizing the results of the paper. 



2 Description of the system 

We study the effect of the difference on restitution coefficients in the segregation phenomenon, 
by means of Molecular Dynamics simulations of a bidimensional granular mixture of two types 
of particles, A and B. Grains are modeled as smooth Inelastic Hard Disks, but differing on 
the normal restitution coefficient that characterizes their inelastic collisions. The two species 
can also differ in their masses tua and ttib, and in their diameters a a and as- The restitution 
coefficient for A- A collisions is ua, for B-B collisions is as. For the interparticle collisions A-B 
we have taken aAB = {(^A + oib)I^- Usually we will consider that B are the most inelastic 
particles (as < (xa)- 

We have taken a fixed total number of particles Nt = Na + Nb , changing the concentration 
of the B particles. For the simulations reported in this paper, we have fixed Nt = 680 disks 
and varied Nb from 10 (that can be considered as a tracer limit) until 160. The disks are 
placed under the action of a gravitational acceleration g pointing downward in a rectangular 
box of width — 50cryi, infinite height, and with the bottom wall oscillating periodically at 
high frequency w and small amplitude A, with a bi-parabolic waveform |23j . Periodic boundary 
conditions are used in the horizontal direction. They avoid the appearance of convective rolls 
by the influence of the walls. Under these conditions, the system reaches a stationary state with 
gradients in the vertical direction |24) . 

Units are chosen such that aA = cr = 1, ruA = 1 and kB = ^- The time scale is fixed by the 
characteristic energy of wall oscillation, that is mA{Aui)'^ — 1. Simulations are performed with 
a fixed value of the gravity acceleration g in order to reduce the parameter space and provide 
a more detailed analysis on the effect of inelasticity. The value of g = 0.15 (in our units), was 
chosen because the associated value of the ratio of the acceleration induced by the vibration 
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versus the gravity, F = {Aw'^)/g can be easily reached experimentally, as it takes the value of 
r ~ 6.67. Finally, the amplitude on the vibration is set to A = O.Olcr. 



3 Inelasticity induced segregation 

In this section we reproduce the main results concerning the macroscopic segregation of the two 
species when they only differ in restitution coefficient, that were already reported in Ref. |22| . 
First, to illustrate the main observed features, we report results of a simulation having a small 
fraction of inelastic particles Nb = 10, and ob ~ 0.7 and the rest nearly elastic: Na = 670, 
aA = 0.98. 

The density profiles of the two species are shown in Fig. [1^, where we plot the number 
density of particles of type A and B: nA{z) and nB^z). The normalization of these quantities 

is such that dz J^'^ dxnA{z) = Na (resp. B). For plotting purposes only, ns is rescaled 
by a factor Na/Nb, so in the case of no segregation both profiles would be identical. Both 
densities have the characteristic shape of vibrofiuidized systems subject to gravity: there is a 
initial density increase due to the abrupt temperature drop caused by dissipation, and at higher 
positions, density decreases again due to gravity [23]. Density exhibits a maximum at z ~ 15a 
where the density n ~ 0.5, so the system cannot be considered as dilute. The density profile 
of the more inelastic particles, B, is plotted as a dotted line in Fig la. Its maximum is located 
at smaller z, indicating that they are closer to the bottom of the container as compared to the 
more elastic ones. The segregation is not complete, though. 

The temperature profiles are also highly inhomogeneous, as shown in Fig. [2^. For both 
species, the temperature presents a initial abrupt drop, but later (after z ~ 20cr) both profiles 
present a linear increase with height. This phenomenon was already observed in one-component 
systems, and it is associated to the energy transport term, — /xVn, related to density gradients, 
that appear in granular fluids [15] . Let us note that the maximum density does not coincide with 
the temperature minimum: there is a shift between these two quantities which is qualitatively 
described by a hydrostatic balance in presence of gravity [25l[T5] . 




Fig. 1. Density profiles for the two species, density of A, ua, (sohd line) and the rescaled density 
of B, ubNa/Nb, (dotted line), (a): both species are inelastic with Na = 670, Nb ~ 10, qa ~ 0.98, 
as — 0.7. (b): A is elastic while B is inelastic and Na ~ 640, Nb = 40, and as ~ 0.5 

The described segregation of species A and B is produced by their different restitution 
coefficients as all the other properties are the same. To study in more detail the effect of the 
difference of inelasticities, we proceed to study the limiting case in which the A particles are 
elastic {aA = 1) and only the B particles are inelastic (and consequently, collisions A-B are 
also inelastic). In this way we also limit the parameter space, allowing to a more detailed 
quantitative study. 

Figures [T)d and[2}D show the density and temperature profiles for such case, where particles 
of type A are elastic {aA = 1) particles B are inelastic (as = 0.5) considering Nb = 40. It 
is observed that the main properties of the profiles are preserved, even the positive slope of 
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z z 

Fig. 2. Temperature profiles for two species, Ta (solid line) and Tb (dotted line), (a) and (b) graphics 
have the same parameters as in Fig. [T] 



Ta despite the A-A collisions are elastic. Partial segregation is again observed, where inelastic 
particles, B, sink to the bottom of the container while elastic ones. A, are majority at upper 
layers of the fluid. 

These results are not surprising in view of the predictions of Ref. |15j . where it is argued 
that the segregation is produced when the particles with different restitution coefficients are 
immersed in a temperature gradient. The gradient is induced by the inelastic collisions, so such 
gradient can be created vibrating a mixture of elastic and inelastic particles. The latter ones 
dissipate some of the energy injected by vibration creating a stationary state. The hydrodynamic 
description of the mixture also contains the dissipative flux — /iVn, and therefore it is expected 
that the hydrodynamic profiles of density and temperature will be equivalent to a full inelastic 
system. 

Figure [3] shows the temperature ratios Tb(z)/Ta{z) for the simulations described in Fig. 
[51 At low densities it was found experimentally [T71IT5] and by employing kinetic theory [TB] 
that such temperature ratio must be constant. However, we find a non constant ratio in the z 
direction that only agrees with the result of [TB] at high z, say for z larger than 30. Inspection of 
Fig. [T] shows that the density at z > 30 is smaller than 0.2, so we can consider that the system 
is in a low density regime. Therefore the approximations of [TB] are valid. In the case where A 
particles are elastic, equivalent predictions for the temperature ratio were given in [26[|27) . 
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Fig. 3. Temperature ratios Tb{z) /Ta{z) for the simulations described in Fig. [21 Simulation results 
(solid circles) and theoretical prediction in j16 (solid line) for Na = 670, Nb = 10, ua ~ 0.98, 
qb = 0.7. Simulation results (open squares) and theoretical prediction in [TS] (dashed line) for A 
elastic, iV^ = 640, Nb = 40, and ob =0.5 
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In order to quantify the segregation, a series of simulations are performed with iV^ ranging 
from 10 to 160, and as between 0.2 and 0.9. Larger values of Nb or smaller restitution coeffi- 
cients lead to clustering as described in ^28^. For each simulation we compute the segregation 
parameter, defined as: 



6 = 1 — J dznA{z)nB{z) j \j \ dzn\{z) J dzn\{z) (1) 

where the nA{z) and nB{z) are the local density, as plotted in Fig.[T] The segregation parameter 
is bounded between and 1. The value 5 = 1 corresponds to complete segregation, as 5 is 1 if 
n^(z) and ub^z) do not overlap. On the contrary, 6 = means complete mixing, as this value 
can only be obtained if nsiz) is proportional to n^(z). 

The results for S are collected in Fig. 0] where the quantity 6 is plotted versus the coefhcient 
as for different values of Nb ■ The fact that S is always non vanishing confirms that the segre- 
gation exists whenever the restitution coefficients are different. Only in the case when as — > 1 
the quantity S approaches 0, limit in which there is no segregation. Note that, for each as, (5 
increases with Nb- The results confirm that segregation is not complete as 6 never gets close 
to 1. In addition, for each simulation, the center of mass of the A and B species are computed, 
Za/B: finding that Za > Zb in all cases. 




The reported inelasticity induced segregation have a microscopic origin as stated in Ref. 
[22] : the most inelastic particles induce locally a region of high density and low temperature, 
resembling a cold droplet that falls in a gravitational field. The droplet is created by the 
dissipation in a way that resembles the clustering instability of the granular gases. Quantitative 
proof of this argument can be found in Ref. |22) , by studying the density and temperature pair 
distribution functions of the A and B particles. 



4 Competition with buoyancy forces 

Up to now, we have considered a mixture where grains of both species have the same masses 
and diameters. The main conclusion so far is that the inelasticity can create an effective force, 
pointing downwards, that sinks the more inelastic particles to the bottom of the container. 
Besides, when densities differ, either because the mass and/or the diameters of the particles are 
different, buoyancy forces appear making lighter grains move to the top. It can be asked whether 
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the inelasticity induced segregation could compensate the buoyancy force experienced by lighter 
B particles. To confirm or deny this idea, a series of simulations with equal sized particles 
{(Tb/o-a — 1), keeping fixed ua, Na, and Nb, but varying mB/iTiA and as is performed. In 
order to determine which force dominates, we compute the center of masses of the A and B 
particles, denoting them by Za, Zb- Their precise definition is: 

^ JPIP^ ^ iTi^, (resp.i.). (2) 

If Za > Zb particles of type A will be, on average, closer to the top of the container than 
particles of type B and viceversa. The ratio of Zb/Za being larger or smaller than one will be 
our criterion for deciding which force actually dominates. 

The results for a a — 1, Na — 640, and Nb — 40 are presented in Fig. \5[ The mass ratio 
mB/iTiA is varied from 0.3 until 1, and the inelasticity of the B particle from 0.4 until 0.9. 
For each value of a^, there is a particular value of the mass ratio, (ms/mA)* at which the 
positions of the two center of masses coincide. The inelastic particles have lower center of mass 
if niB/rnA > (jnBlmA)* , and therefore sinking due to inelasticity wins to the buoyancy force 
in this range. On the contrary, buoyancy force dominates if ms/m^ < (ms/myi)* and lighter 
and inelastic particles go up. 

Figure also presents the value of the critical ratio {mB/rnA)* as a function of cub'- more 
inelastic particles have a higher tendency to sink and therefore a stronger mass ratio is needed to 
compensate the inelasticity induced segregation. Note that when as — ^ 1, the positions of both 
center of masses diverge because a vibrated elastic medium does not have a stationary state. 
There is an indication of such divergence in the steep slope of the dashed line in the upper 
left panel of Fig. [5] The segregation parameter 5 does not vanish for any value of mB/mA, 
indicating that there is no complete mixing even at the value of mB/mA = (mB/mA)* , where 
both center of masses coincide. The value of S (not shown here), however, is minimum at this 
precise mass ratio. Also, S does not reach 1 either, so no complete segregation is observed. 

A second series of simulations were done with aA = 0.9 and are presented in Fig. [S] The 
qualitative features are similar to the the previous case (a a = 1-0), except for some remarks. 
The first one in the center of mass of species A is about constant, independent of the mass 
ratio or the dissipation parameter of particles B. This is so because the particles A are inelastic 
and dissipate a great deal of energy, forming a very dense system. This density is near the 
close packing, as it will be illustrated in the forthcoming figures. The second remark is that, 
as opposite to the a a = 1 case, the center of masses do not diverge when as — > 1 due to the 
finite inelasticity of the A particles. Finally, as a a = 0.9, the ratio ob/cka can exceed one. In 
particular, we can make the B particles elastic by choosing ob = 1/0.9 = 1.1111.... We have 
done so and included the results in the bottom part of Fig. [5] As expected, particles B must 
be heavier than A in order to compensate the effective upwards force produced because they 
are more elastic. In this particular case the value of (mB/mA)* — 1.17. 



5 Competition with the Brazil nut effect 

When grains are mixed such that they differ in size, a remarkable phenomenon takes place: the 
larger grains rise to the top even if they are denser. This phenomenon, called the Brazil nut 
effect (BNE), has been widely studied in order to identify the different mechanisms that produce 
it and also to characterize the range of parameters where it takes place. As mentioned in the 
introduction, it has been found that besides the normal BNE, it is also possible that, under 
certain conditions, the large particles sink in the so-called reverse Brazil nut effect (RBNE). 
In its simpler form, the BNE is studied in the neutrally buoyant case, where the mass ratio 
between the species equals their volume ratio. 

The results of the previous sections suggest that if the larger particle is also more inelastic 
the inelasticity induced segregation can compete with the BNE. We have performed a series of 
simulations of a neutrally buoyant mixture (i.e., ms/m^ = (cb/cta)^) where we have varied the 
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Fig. 5. Left: Position of the center of mass of A (open squares) and B particles (solid circles) as a 
function of the relative mass of B particles mB/mA for different values of the ratio of the restitution 
coefficients, as /a a (indicated in the plot). The locus of the mass ratio at which the two centers of mass 
coincide is indicated by a thick dashed line. Right: Critical mass ratio (ms/mA)* at which the two 
centers of mass coincide as a function of the ratio of the restitution coefficients as/oA- The regions in 
the parameter space where B particles move to the top or to the bottom are indicated by the respective 
labels in the plot. Simulation parameters are fixed to Na ~ 640, Nb = 40, (Jb/(^a = 1, and qa = 1 
(top) and a A ~ 0.9 (bottom). 



size of B particles (Tb and the restitution coefficients a a and as, while the number of particles 
are fixed to Na — 640 and Nb = 40. This parameter space is huge as we have restricted to B 
particles that are not too large; specifically we have considered 1 < crs < 2.5. 

First, we consider the usual BNE case with aA = as and as = 1.5, 2.0 (that is, the two 
species only differ in their sizes). At low dissipation {aA = as = 0.95) there is no BNE but 
rather the two species mix with B particles having a narrower distribution near the center 
(see Fig. [S]). When the dissipation is higher {ua = ctB = 0.9,0.85) the larger particles sit on 
the top as in the usual BNE (see Fig. The segregation is partial when as — 1-5 but it is 
complete when as = 2.0, differently as in the other cases, where segregation was always partial. 
In summary, we have checked that the BNE can be obtained with moderate size difference as 
long as the inelasticity is strong enough. 
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When the two inelasticities differ (with B particles being larger in size and more inelastic), 
it is expected that the BNE could be compensated by the inelasticity induced segregation. A 
series of simulations with a a = 1, ctb = 2.0, and different values of as show that the inelasticity 
induced segregation is the dominant mechanism and B particles sink (that is, RBNE is found) 
for large inelasticities and when the elastic limit is approached both species mix (see Fig. [T]) . 
However, the elastic limit could not be achieved because, when both species are elastic, the 
vibrated granular gas has no stationary state. 

A different situation is found when aA = 0.9 and gb = 1-5, where a transition from RBNE 
to BNE is obtained with coexistence of both states. The density profiles for different values of 
as are presented in Fig. [5] When both inelasticities are equal or the B particles are slightly 
more inelastic than A particles, a clear BNE is found. But when the B inelasticity is increased, 
an important fraction of B particles sink showing a coexistence of BNE and RBNE. In fact, 
a "sandwich" configuration is obtained with B particles being in the bottom and top layers 
while A particles are in the middle region. The segregation, as is common in all the examples 
presented here, is not complete and particles have a partial mixing. An stronger separation has 
been found experimentally |29j in a granular mixture of glass and bronze particles vertically 
shaked. These authors found, however, complete segregation. 

If we define the transition from BNE and RBNE as the dissipation when Za = we find 
a value for as — 0.71. In other words, when gb = 1-5, a dissipation as < 0.71 makes the B 
particles to sink. As mentioned before, the parameter 5, defined in Eq. ([T]), takes its minimum 
value: 5 = 0.053. 

The coexistence of BNE and RBNE could be a transient configuration due, in part, to the 
high compacity of the small grains in the middle region that limits the mobility of the large 
particles. They could not migrate between the bottom layers and the top ones. A dynamic 




Fig. 6. Left: Top: Density profiles of A particles (solid line) and B particles (dashed line) in a neutrally 
buoyant mixture with an = 2.0, Na ~ 640, Nb = 40, and oa = as = 0.95 (left) and oa = as = 0.85 
(right). Bottom: Configuration of the grains in the two systems. 
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Fig. 7. Density profiles of A particles (solid line) and B particles (dashed line) in a neutrally buoyant 
mixture with qa = 1, (J_b = 2.0, Na ~ 640, Nb ~ 40, and different values of as (indicated in the 
figures) . 



analysis, that is beyond the purpose of this article, could lead some insight into the coexistence 
of both effects. 

An inspection of the particle configurations with as = 0.3 where both RBNE and BNE 
coexist, shows that B particles develop a tendency to form pairs and small clusters (see FiglH]). 
This phenomenon was already reported in the case of particles of equal size and it was quan- 
titatively characterized in terms of the pair distribution functions j22j . There is not only a 
macroscopic segregation described in terms of the density profiles and center of mass positions, 
but also a microscopic segregation with a tendency of the more inelastic particles to aggregate 
in small clusters. In this case, the pair distribution function of the AA particles resembles the 
structure of a hexagonal crystaline structure, while the BB pair distribution function has a 
very high peak at ctb, followed by a much smaller peak and hardly any structure for distances 
larger than 2a b- 

We have shown that inelasticity induced segregation can compete with both buoyant forces 
and the BNE, separately. If a mixture of grains is considered such that B particles have the same 
mass as A particles, but have larger size, both buoyancy forces and the BNE push B particles 
to the top. The simulation results presented in Fig. [TUl show that for ctb > 1.9 buoyancy and 
BNE cannot be compensated by inelasticity and the large grains go always to the top. When 
as ^ 1-9 inelasticity induced segregation can balance the the combined effect of buoyancy and 
BNE and higher inelasticities are needed to compensate their effects for larger B particles. 

The temperature profiles present the same qualitative behavior as in the case of a mixture 
where grains differ only in their restitution coefficient (see Fig. 1111) . That is, an initial abrupt 
drop, followed by a slight linear increase with height. The temperature ratio Tb{z)/Ta{z) shows 
that B particles are always colder than A particles but it does not agree with kinetic theory 
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Fig. 8. Density profiles of A particles (solid line) and B particles (dashed line) in a neutrally buoyant 
mixture with ua = 0.9, ob = 1-5, Na = 640, Nb ~ 40, and different values of as (indicated in the 
figures) . 



predictions for dilute and moderately dense granular gases [16) as we work in a high density 
regime, close to close packing in some of our simulations. 



6 Conclusions 

The main conclusion of the present paper is that different restitution coefficients alone create 
segregation in a binary mixture vertically vibrated. This segregation mechanism can compete 
or enhance other segregation mechanisms as buoyancy or the Brazil nut effect. The restitution 
coefficients must be considered, in addition to the usual material properties (mass ratio and 
diameter ratio) , in order to describe accurately the segregation process in a granular mixture. 
The effect of the inelasticity is such that the most inelastic particles sink to the bottom of the 
container while the less inelastic ones raise to the top. Segregation is not complete, however, 
but only partial. The density profiles of each species show a maximum which are located in 
different positions depending on the inelasticities, but they have an important overlap. 

The segregation effects presented here also appears by vibrating a mixture of elastic and 
inelastic particles. Again inelastic particles migrate to the bottom of the container and elastic 
one prefer the upper part. This phenomenon allows to make a systematic study of the mechanism 
reducing the parameter space. 

If, besides the difference of inelasticities, the two species differ in mass density buoyancy 
forces enter into play. If the inelastic particles are less dense, the two mechanism compete 
and a transition line in the parameter space is found that distinguish the cases where the 
inelastic particles have lower or higher center center of mass compared to the elastic ones. 
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Fig. 9. Configuration of the grains with ua = 0.9, as = 0.3, as = 1.5, Na = 640, and Nb = 40. The 
large and more inelastic grains aggregate in small clusters. 




Fig. 10. Left: Position of the center of mass of A (open squares) and B particles (solid circles) as a 
function of the relative size of B particles as/cTA for different values of as (indicated in the plot). 

The locus of the mass ratio at which the two centers of mass coincide is indicated by a thick dashed 
line. Right: Critical size ratio ((Jb/cta)* at which the two centers of mass coincide as a function of the 
inelasticity of B particles, qb/qa. The regions in the parameter space where B particles move to the 
top or to the bottom are indicated by the respective labels in the plot. Simulation parameters are fixed 
to OA = 1.0, Na = 640, Nb = 40, and tub = 1. 



More inelastic particles need higher buoyancy forces to compensate the inelasticity induced 

segregation mechanism. 

When the two species have the same density (that is, tlicy are neutraUy buoyant) but differ 
in their size, the Brazil nut effect appears: large particles have the tendency to go to the top. 
When the large particles are also more inelastic, the two mechanism compete. The simulations 
show that there is a range of parameters where a large number of the large and inelastic grains 
arc simultaneously at the top and at the bottom: Brazil nut effect and the reverse Brazil nut 
effect coexist. The relative proportion of the particles in the two regions depend on the control 
parameter. 
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Concerning the temperatures, most dissipative particles have a lower temperature than the 
most elastic ones, both having the characteristic shape of vibrofluidizcd system (fast cooling 
away from the moving boundary followed by a heating that grows linearly with the distance). 
Their rate, however, does not agree with kinetic theory predictions for dilute and moderately 
dense granular gases. 
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